Phase diagram of the one-dimensional Hubbard model 
with next-nearest-neighbor hopping 
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We study the one-dimensional Hubbard model with nearest-neighbor and next-nearest-neighbor 
hopping integrals by using the density-matrix renormalization group (DMRG) method and Hartree- 
Fock approximation. Based on the calculated results for the spin gap, total-spin quantum number, 
and Tomonaga-Luttinger-liquid parameter, we determine the ground-state phase diagram of the 
model in the entire filling and wide parameter region. We show that, in contrast to the weak- 
coupling regime where a spin-gapped liquid phase is predicted in the region with four Fermi points, 
the spin gap vanishes in a substantial region in the strong-coupling regime. It is remarkable that a 
large variety of phases, such as the paramagnetic metallic phase, spin-gapped liquid phase, singlet 
and triplet superconducting phases, and fully polarized ferromagnetic phase, appear in such a simple 
model in the strong-coupling regime. 
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I. INTRODUCTION 

For several decades, quasi-one-dimensional (ID) mate- 
rials have been one of the major subjects of research in 
the field of condensed matter physics^! 3 - A standard de- 
scription of such materials is the ID Hubbard model.— 
The simplest case with the cosine dispersion (nearest- 
neighbor hopping only) was solved exactly by Lieb and 
Wu via the Bethe ansatz^ The low-lying excitations were 
also understood well as the Tomonaga-Luttinger liquid 
(TLL),— where the renormalization group technique and 
bosonization method have been usedi* 5 - However, modi- 
fications of the ID Hubbard model are often required for 
realistic descriptions of the materials. In general, such 
modifications (even if they are small) make the analy- 
ses much more difficult since the correlation effects are 
strong in low-dimensional systems. Thus, even in the ID 
systems, our knowledge is still far from being complete. 

One of the typical modifications is to add a next- 
nearest-neighbor hopping term in the Hamiltonian, which 
brings a sort of frustration to the spin degrees of free- 
dom of the system as well as some coupling between spin 
and charge degrees of freedom. In the past, this model 
has been extensively studied and some distinctive fea- 
tures, which are absent in the simple ID Hubbard model, 
have been found. At half filling, the system has three 
phases: one is a Mott insulating phase with 2fcp spin- 
density-wave (SDW) correlation (which occurs when the 
spin frustration is small); the others are a spin-gapped 
insulating phase with incommensurate spiral correlation 
and a spin-gapped metallic phase for sufficiently large 
spin frustration . 9 ' 10 i 11 i 12 Away from half filling, the ex- 
istence of ferromagnetism has been shown analytically 
in some limiting cases for infinite strength of the cou- 
pling ) 34 ' 35 ' 36 which has been confirmed numerically for 
finite but large enough strength of the coupling i 13 ' 14 ' 15 
Also, it has been pointed out that, although a weak- 



coupling analysis leads to only a spin-singlet supercon- 
ducting phase with finite spin gap^ previous density- 
matrix renormalization group (DMRG) studies suggest 
that the spin gap vanishes for large enough coupling 
strengths when the next-nearest-neighbor hopping is pos- 
itive and large . 10 ' 18 Moreover, a spin-triplet supercon- 
ducting phase has been shown to exist at quarter fill- 
ing* 1 ^ As just described, the present system has many 
phases unparalleled in other ID strongly-correlated elec- 
tron systems; i.e., our modified ID model can involve a 
variety of physical phenomena. In particular, to detect 
their phases is of particular interest in the light of recent 
proposals to realize a Hubbard model of fermions on an 
optical lattice.— 

There are some relevant materials to the ID Hub- 
bard model with next-nearest-neighbor hopping. One is 
the quasi-lD organic conductor (TMTSF) 2 X [X=PF 6 , 
CIO4], the so-called Bechgaard salt. 2 i*22 This material 
exhibits a rich phase diagram upon variation of the pres- 
sure and temperature. At low temperatures, the phase 
changes in the order as the spin-Peierls insulator, an- 
tiferromagnetic insulator, spin-density-wave (SDW) in- 
sulator, superconductor, and paramagnetic metal, with 
increasing pressure. So far, experimental evidence that 
the superconducting state is in the spin-triplet channel 
has been piled up i 23 ' 24 Theoretically, it has been pro- 
posed that a triangular lattice formed by the hopping 
integrals makes the ferromagnetic ring-exchange mech- 
anism relevant, which in turn leads to the spin-triplet 
superconductivity. 25 

Another relating system is a newly synthesized copper- 
oxide compound P^Ba^uyOis-^.— This material con- 
sists of both the single CuO chains (as in PrBa2Cua07) 
and double CuO chains (as in PrBa2Cu40s), and those 
chains are separated by insulating Cu02 plains. It has 
been reported that the double chains turn into a super- 
conducting state below T c ~ 10 Ki2! So far, some nu- 
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merical studies have been carried out; on the basis of the 
TLL theory, a weak-coupling phase diagram has been 
obtained in the d-p double chain model . 28 > 29 Also, in 
a reduced single-band double chain model, it has been 
suggested that the superconducting gap has an extended 
s-wave-like form, which does not contradict with the ex- 
perimental results^ Relevance of charge fluctuations has 
also been discussed i 31 i 32 

In this paper, we study the ID Hubbard model with the 
next-nearest-neighbor hopping. We calculate the TLL 
parameters, spin gap, and total-spin quantum number 
by using the DMRG method and Hartree-Fock (HF) ap- 
proximation. Based on the results, a detailed phase dia- 
gram as a function of band filling and hopping integrals 
is determined in both weak-coupling and strong-coupling 
regimes. Surprisingly, the phase diagram of the model 
contains a large number of distinct phases in the strong- 
coupling regime although the model is quite simple. We 
hope that the present investigation will contribute to bet- 
ter understanding of the ID strongly-correlated electron 
systems. 

This paper is organized as follows. In Sec. II, we define 
the ID Hubbard model with the next-nearest-neighbor 
hopping and study the model in the noninteracting case. 
In Sec. Ill, we discuss how the Hartree-Fock approxima- 
tion and DMRG method are used to calculate the TLL 
parameter. In Sec. IV, we present the calculated results 
and obtain the phase diagrams of the model based on 
the numerical results. Section V contains summary and 
conclusions. 



II. MODEL 

We consider the ID Hubbard model with the next- 
nearest-neighbor hopping, which is defined by the Hamil- 
tonian 

H = ^I]( C l + la C -+ H - C -)+^^(ci +2CT C^+H.C.) 

i,cr i,cr 

+ U ^n^riii, (1) 

i 

where cJ CT (cj CT ) is the creation (annihilation) operator of 
an electron with spin a at site i, and n; L(J = C ig Ci a is 
the number operator. t\(> 0) and t 2 are the nearest- 
neighbor and next-nearest-neighbor hopping integrals, 
respectively, and U is the on-site Coulomb interaction 
[see Fig. [T{a)]. 
The dispersion relation is given by 

Bk = 2ii cos ka + 2t 2 cos 2fcet, (2) 

where k is the wave number and a is the lattice constant; 
we set a = 1 hereafter. The Ijare band width is esti- 
mated as W = 2<i + 4|i 2 | + for \t 2 /h\ > 1/2 and 
as W = 4ii for 1 1 2 /*x | < 1/2. The ground-state phase 
diagram in the noninteracting case (U = 0) is shown in 
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FIG. 1: (a) Schematic representation of the lattice structure 
of the ID Hubbard model with the next-nearest-neighbor hop- 
ping and (b) its U = phase diagram, where the thick solid 
curves separate the two regimes. The inset shows the qual- 
itative behavior of the band dispersion et, where the Fermi 
level is indicated by the thin line. 

Fig-Gib). For \t 2 /h\ < cos 2 [(2-n)7r/2]/ sin 2 [(2-n)7r] (n 
is the band filling) , the system has two Fermi points and 
the physical properties at low energies are qualitatively 
the same as a system with t 2 = 0. On the other hand, for 
|*2 /*i | > cos 2 [(2 — ti)tt/2}/ sin 2 [(2 — n)7r], there are two 
branches, namely four Fermi momenta ±fcpi and ±&F2 
(\k F2 \ > \k F1 \). In this discussed in Sec. IV. A, 

the Fermi surface can be mapped to that of a two-leg 
Hubbard ladder model at weak coupling. We designate 
the critical boundary at which the Fermi surface splits 
into four points as the Fermi-point (FP) boundary and 
the FP boundary is characterized by a condition kp\ — 
(or &p 2 = tt)- Hence, the model ([T]) has to be dealt with 
as a two band system within the TLL theory. Note that 
the parameter region n > 1, t 2 > is exactly equivalent 
to the region n < 1, t 2 < under the particle- hole trans- 
formation. By the same token, the region n < 1, t 2 > 
equals to the region n > 1, t 2 < 0. We therefore consider 
only the region < n < 1 for both positive and negative 
values of t 2 . 



III. METHOD 

The low-energy properties of TLL are characterized by 
a few quantities; most notably, the TLL parameter K p 
determines the long-range behavior of various correlation 
functions in the metallic TLL ground state. It has how- 
ever been recognized that the numerical calculation of 
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K p for an arbitrary strength of correlations is very diffi- 
cult. Recently, one of the authors has succeeded to over- 
come this difficulty,— where a simple and stable method 
for calculating K p with the DMRG technique in single- 
band ID systems is proposed. In this section, we extend 
the method to the two-band systems and check the per- 
formance of the method by comparing the results with 
those obtained by the HF approximation which is known 
to provide a good estimation of K p in the weak-coupling 
regime. 



A. HF approximation 

It is known that the small-C/ perturbative estimation of 
K p is feasible for U < W/2 in the ID Hubbard models 
One of the authors applied this perturbative method to 
a two-leg Hubbard ladder model with four Fermi points 
and confirmed that it gives quantitatively reliable results 
in the weak-coupling regime (U < W/A)2^ Since the low- 
energy physics of the two-leg ladders is equivalent to that 
of our model defined by Eq. |l])r^ we may naturally ex- 
pect the perturbative estimation to be applicable to our 
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In the TLL theory, the critical exponent K p is given 



1 



where x is the charge susceptibility defined as 

_ a = I d 2 E gs {n) 
L dn 2 ' 

and D is the Drude weight defined as 



D 



7T d 2 E gs {0) 
L dd> 2 ' 



(3) 



(4) 



(5) 



Here, L is the number of lattice site, n is the band filling, 
£" gs is the ground-state energy, and <\> is the magnetic 
fluxi^ Within the first-order perturbation expansion, E gs 
can be determined as 



Egs = E 



UL 



where Eq is the ground-state energy of the corresponding 
noninteracting system. We then obtain 
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(7) 



where xo and Do are the charge susceptibility and Drude 
weight of the noninteracting system, respectively. A sim- 
ple expression for K p is therefore obtained as 



Kn 



Uxo 



(8) 



Note that this scheme is equivalent to the HF approxi- 
mation^ 
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FIG. 2: Values of K P (L) calculated by the DMRG method 
and plotted as a function of the inverse system size 1/L. Solid 
lines are the polynomial fits to the data for finite-size scaling 
analysis. Diamonds indicate the results calculated from the 
Hartree-Fock approximation. The upper and lower panels 
show the result at fo/ti = —3 and tz/ti = 3, respectively. 
U/ti = 2 and n = 0.2 are assumed in both panels. 



B. DMRG 

With the DMRG method, the TLL parameter K p is 
most generally obtained from the long-range decay of the 
density-density correlation. The density-density correla- 
tion function is defined by the ground-state expectation 
value 



1 L 

C m (r) = -J2((ni + rni)-{n l+r ){ ni )). 



(9) 



1=1 



When the system has two Fermi points, it is known that 
the asymptotic behavior is given by 



C NN (r) 



K p Acos(2k F r) _ 3/2 
In ' (r 



(io) 



(6) where A is a constant . 37 ! 41 We can extract K p via the 



Fourier transformation of Eq. 



C NN (g) = A^ e -^C NN (r), 
i=i 



(11) 



where < q < 2n. From the derivative at q = 0, one 
finds the expression 



K p {L) = |G NN (I 



K p = Urn Kp(L). (12) 

L — >oo 



for the thermodynamic limit. It has been demonstrated 
that the value of K p can be determined quite accurately 
by using Eq. (fT2")l with the DMRG method for the single- 
band Hubbard model. Thus, for the precise estimation, 
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we need to calculate the density-density correlation func- 
tion directly in the Fourier space; see Ref. [53 for further 
details. 

Let us now apply this scheme of estimating K p to a 
system with four Fermi points. We then have to assume 
the asymptotic behavior of the density-density correla- 
tion function. Here we assume the behavior 



C NN (r) 



2K p Bcos[2(fc F1 - k F2 )r] 



, (13) 



in analogy with the case of two coupled chains,— because 
the calculated low-energy excitation spectra of our model 
are similar to those of the two coupled chains.^ We thus 
obtain 



K P (L) = ^ (| 



K p = lim K p (L), (14) 

L — >oo 



as a substitute for Eq. (jl2|) . In principle, one may cal- 
culate Eq. ([H]) in the same way as Eq. (JT3J) . However, 
the discarded weight in the DMRG calculation increases 
rapidly with increasing I^Ail, so that the calculation 
must be carried out more carefully. 

In this paper, we apply the open-end boundary con- 
ditions for precise DMRG calculations^ We keep up to 
m ~ 4500 density-matrix eigenstates in the DMRG pro- 
cedure and extrapolate the calculated quantities to the 
limit m —* oo. We also use several chains with lengths 
L = 40 to 240 and then perform the finite-size scaling 
analysis based on the size-dependence of the quantities. 
In this way, we can obtain the quite accurate ground state 
with an accuracy of AE gs /L ~ - lO -5 ^. In Fig. H 
we demonstrate the finite-size scaling analysis for (a) the 
two- and (b) four-Fermi point cases. For both cases, one 
can see the systematic extrapolation of K p to the ther- 
modynamic limit L — > oo. We also find that, at least in 
this coupling strength U/ti — 2, a good agreement is ob- 
tained between the extrapolated values of K p obtained 
from the DMRG data and the corresponding values of 
K p obtained from the HF approximation. 
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FIG. 3: Phase diagram of our model in the weak-coupling 
limit U = + . We follow the notation of the symbol CaS/3 of 
Ref.-£. 



Daul and Noack adapted the weak-coupling analysis of 
Ref. for the analysis of the Hamiltonian ([T])^ 

In Fig. [H we show the phase diagram for ^Ail < 3 
and n < 1 in the weak-coupling limit U = + . A notation 
CaS/3 denotes a phase with a gapless charge modes and 
[3 gapless spin modes, where a and (3 are integer values 
from to 2. Generally speaking, a metallic phase with 
four (two) Fermi points is characterized by CISO (ClSl). 
Also, a spin-gapped liquid phase CISO appears around 
the FP critical boundary due to the van Hove singularity 
of the model. Note that the TLL parameter has the value 
K p = 1 in nearly all the metallic regimes, except at the 
FP critical boundary and on the C0S1 line, where we 
have K p = 1/2. 



B. Small U 

Let us turn to the small-{7 perturbative regime where 
we choose the coupling strength U — 2t\. In this regime, 
the HF estimation of K p [Eq.©] is expected to give a 
good approximation. When the system has two Fermi 
points, we can simply obtain 



IV. RESULTS 

A. Weak-coupling limit 

Let us first consider the phase diagram within a weak- 
coupling analysis. Balents and Fisher have obtained the 
weak-coupling phase diagram of the two-leg Hubbard lad- 
der model using the renormalization group technique and 
bosonization method. A - Their analysis can be universally 
applied to a system with four Fermi points. In addition, 
Fabrizio pointed oufeil that the low-energy physics of the 
two-leg Hubbard ladder model is exactly the same as 
that of our model |T]) via a simple mapping of the Fermi 
points, i.e., ±fcpi — > =Ffcp and ±k?2 — * T^f wnere 
(±fc F ) are the Fermi points for the antibonding (bond- 
ing) band of the two-leg Hubbard ladder model. Then, 



which leads to 



TTVp 



TTVp + U 

where the Fermi velocity is 

v-p = 2t\ sin(7T7i/2) — 4^2 sin(7rn). 



(15) 



(16) 



(17) 



This expression (flU)) is universal to the TLL with two 
Fermi points. On the other hand, when the system has 
four Fermi points, we obtain 



-1 7T \Vf\VF2\ 



2 |ufx| + |w F2 | 



(18) 



5 



(a) 



" 0.9^ / ]/IJ\ 
J 0.9/ /III > 


i 




1 K p =Vg5 ■ 

I \0.9 


*""""* 0.95/ / 


i 


\\\p-8 \^ ■ 



-4 -2 2 4 

t 2 /ti 




-4-3-2-101234 

*a/*i 



FIG. 4: Contour map of the TLL parameter K p in the 
n—til%\ plane for the weak-coupling interaction U = 2t\. The 
results are obtained from (a) the HF and (b) DMRG calcula- 
tions. The thick line represents the FP critical boundary. 



When the system has two Fermi points, the ground 
state can basically be presumed to be a standard ID 
TLL. In both contour maps, K p becomes larger as *2 
decreases for fixed n. This is because the effective in- 
teraction parameter U/vf is reduced with decreasing ti- 
In particular, K p varies drastically around t% ~ except 
when n ^ 0.1 and n > 0.9, which is due to the rapid 
change of the inverse Fermi velocity Wp 1 ; for example, 
the Fermi velocity is estimated as vf = \plt\ — 4*2 at 
quarter filling (n — 1/2). Note that K p = 1/2 is reached 
in the limits n — > and 1 and that K p > 1/2 everywhere 
else. 

As soon as the Fermi surface splits from two into four 
points, K p drops (almost) discontinuously to 1/2. When 
the system has four Fermi points, v F goes to zero in 
the limits vfi — > (&fi — * 0) or uf2 — * (&F2 — > 7r) 
[corresponding to the diverging density of states] . Con- 
sequently, the effective interaction parameter diverges, 
U/v F — > co, and the strong-coupling value K p = 1/2 
is produced. This behavior is similar to that of the ID 
single-band Hubbard model in the limit of n — > 0. Then, 
K p increases rapidly as the parameters are away from 
the FP boundary line and gets closer to 1 in the limit of 
| £2 1 — ► 00. In the two-band model, the criterion for the 
dominant superconducting correlation is K p > 1/2. We 
thus find that the superconducting correlation is the most 
dominant in the entire region with four Fermi points. 

C. Large U 



after some calculations, where i>fi (i>F2) is the Fermi ve- 
locity at the momentum &fi (&F2) given as 

«F(l,2) = 2ti sin(/c F (i,2)a) - 4i 2 sin(2fc p(1>2 )a). (19) 
Substituting Eq. (JTSJ) for Eq. ©, we obtain 



U 



with using an effective Fermi velocity 



^F 



«FlfF2 



1>F1 



\VF2\ 



(20) 



(21) 



In Fig. 31 we show the contour maps of K p at U = 
2*i in the parameter space of * 2 /*i and n, which are 
calculated with (a) the HF approximation and (b) the 
DMRG method. The thick line represents the FP critical 
boundary and the thin lines form the contour map of K p . 
We find that the quantitative agreement between the two 
phase diagrams is pretty good, which means that the HF 
scheme is still valid for this interaction strength U/ti = 
2. In the entire region of the phase diagram, we find 
K p > 1/2 and thus the ground state may be described as 
the TLL. However, the deviation of the HF data from the 
DMRG data is relatively large around *2 ~ and n»l, 
where the umklapp scattering becomes dominant. 



Let us now consider how the small-C/ contour map 
is affected by increasing the strength of the on-site 
Coulomb interaction. For large enough coupling U, it 
has been found that there is an extensive ferromagnetic 
phas o 13 i 14 ' 15 and that the spin-triplet superconducting 
state is stabilized at the intermediate filling n ~ 0.5 
when the next-nearest-neighbor hopping is large enough, 
*2/*i ^ 2»i£ A breakdown of the TLL state was also re- 
ported for the ID Hubbard and t— J models with the 
next-nearest-neighbor hopping integrals ; 46 ' 47 where the 
latter model is essentially the same as our model ([1]) in 
the strong-coupling regime. 

We here study the case of U — 10*i as a typical in- 
teraction strength for realistic strongly-correlated elec- 
tron systems. In Fig. we show the contour map of K p 
obtained by (a) the HF approximation and (b) DMRG 
method. We find that, at low densities (n < 0.4), the 
agreement between the two contour maps is qualitatively 
good, while at intermediate to high densities (n > 0.4), 
the situation seems to be totally different. In fact, the 
HF scheme is no longer appropriate and therefore the spin 
and charge fluctuations have to be taken seriously into 
account beyond the usual weak-coupling picture. We will 
thus proceed to a discussion based on the DMRG contour 
map in the following. 

We first note that, as far as the system has two Fermi 
points, the basic properties are qualitatively the same as 
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FIG. 5: Contour map of the TLL parameter K p in the 
n—til%\ plane for the strong-coupling interaction U = 10ti. 
The results are obtained from (a) the HF and (b) DMRG cal- 
culations. The thick line represents the FP critical boundary. 



those of the weak-coupling regime. Thus, the behavior 
of K p at U = lQt i is still similar to that at U = 2t\ 
although the value of K p becomes relatively small. We 
also note that the FP boundary is (slightly) shifted to- 
ward the smaller ^1 direction due to renormalization of 
the band structure at U > 0, as was pointed out in Ref. 
[HI At the same time, the FP boundary line is some- 
what blurred because of some strong quantum effects; 
the change in K p at the FP boundary is still sharp but 
no longer discontinuous as in the weak-coupling limit. 

Let us turn to the case with four Fermi points. The 
ground state is affected drastically by the (strong) inter- 
action strength. Unlike in the small-C/ contour map, we 
find that there is a substantial region with K p < 0.5 
around half filling [denoted by the shadowed area in 
Fig. [5(b)]. According to the TLL theory, this value of 
K p is possible only when long-range repulsive interac- 
tions are included in the model. Therefore, the ground 
state in the shadowed regime would no longer belong to 
the general class of the TLLs. This is consistent with the 
previous DMRG results . 16 ! 47 As discussed in the next 
subsection, this non-TLL-like regime consists of a spin- 
gapped phase and a paramagnetic phase with strong fer- 
romagnetic fluctuations. Also, it is interesting to note 
that the TLL parameter remains constant K p ~ 0.5 in a 
wide region of the phase diagram near the FP boundary 
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FIG. 6: Calculated values of the spin gap (solid symbols, left 
axis) and TLL spin exponent (open symbols, right axis) at (a) 
tz/ti = -1, (b) 1, (c) -2, and (d) 2. U/ti = 10 is assumed. 
The solid line denotes K a = 1 and the vertical dotted lines 
indicate the critical boundaries between the presence and ab- 
sence of the spin gap. 



at low densities. 

Furthermore, it is particularly worth noting that K p 
seems to be enhanced significantly at the FP boundary 
near half filling. [In this area, precise evaluation of K p in 
the thermodynamic limit l/L — ► is rather hard because 
\dK p (L)/d(l/L)\ increases with decreasing the inverse 
system size l/L.] If the definition ([3]) could be still valid 
in this region, K p become quite large > 1: for the van 
Hove singularity, the charge susceptibility \ m Eq. ((4]) 
can diverge and the Drude weight D in Eq. ([5]) must 
remain finite. This may be associated with the C1S0 
phase attributed to the van Hove singularity at kp\ = 
(or Uf2 = 7r) in the weak-coupling phase diagram. This 
is consistent with a prediction that the superconducting 
fluctuations increase with increasing the difference be- 
tween |upi| and |t>F2|, as was suggested in Ref. 



D. Spin gap 

For more elaborate studies of the region with four 
Fermi points, we consider the spin degrees of freedom 
in the strong-coupling regime U — 10ii . Of particular 
interest here is the presence or absence of a finite en- 
ergy gap in the spin excitation spectrum. We thus evalu- 
ate the spin gap defined by an energy difference between 
the first spin-triplet excited state and the singlet ground 
state: i.e., 



lim [E m [N,\)-E m [N,ti)l 



(22) 



where E gs (N,S z ) is the ground-state energy for a given 
number of electrons N and z-component of the total spin 
S z . It is however known that, for some parameter values, 
the spin gap becomes too small to figure out if it remains 
finite, e.g., A s < 10 -3 £i. For verifying the presence or 
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absence of the spin gap, we then calculate the TLL spin 
exponents, which is given by 

L^oo Z * — ' 
kl 

where Sf = n,f — n^. We should find that the spin expo- 
nent takes the value K a — in the spin-gapped phase and 
K a = 1 everywhere else in the thermodynamic limits 
However, for finite-size systems, the situation is not so 
simple. In practice, in the spin-gapless phase, one cannot 
expect to find the value K a — ► 1 due to the logarithmic 
corrections. It is known that the logarithmic corrections 
vanish at which the spin gap opens, in analogy with the 
dimerization transition in the J\-J% model [see Eq. (|24[) 
below] j*i Also, in the spin-gapped phase, if the spin gap 
is very small, the convergence of K a to will obviously 
occur only for very large systems. Therefore, we here 
determine the critical point at which the spin gap opens 
by adopting the condition that the value of K p crosses 
1. This method was first proven to be useful in Ref.— In 
Fig. [SJa)-(d), the spin gap and TLL spin exponent cal- 
culated by the DMRG method are plotted as a function 
of the band filling n for various t2/t\ values. 

For U ^> ti, t2, our model (JTJ) at half filling can be 
mapped onto a Heisenberg model 

i i 

with Ji = At\/U and J x = U%/U. This model has been 
extensively studied both analy ticall y 49 i 50 ' 5 1 1 52 and nu- 
mericallyi 10 ' 18 It has been found that the spin gap opens 
when J2/J1 > 0.241; the ground state is of a dimerized 
zigzag-bond state for 0.241 < J2/J1 < 0.5 and of the 
Majumdar- Ghosh state with incommensurate spiral cor- 
relations for J2/J1 ^ 0.5. For example, the spin gap was 
estimated to be A s ~ 0.25Ji at J2/J1 ~ (ia/ii) 2 = 1 in 
the previous DMRG study^ This value is comparable to 
our estimation A s /ii ~ 0.05 for £2/^1 = 1 and U = 10t\. 
Also, the spin gap is of an exponential dependence on 
J2/J1 as A s oc exp ( — const, x J2/J1) for large J2/J1 
values, which is consistent with a very small spin gap 
A s /h ~ 0.0005 for t 2 /h = 2 and U = 10ii obtained in 
our calculations. 

Let us turn to the evolution of the spin gap upon 
doping. In the weak-coupling phase diagram, a metallic 
phase with four (two) Fermi points is simply character- 
ized as the spin-gapped (gapless) TLL. However, for large 
U values, it is difficult to speculate the n-dependence of 
the spin gap because of the competition between the anti- 
ferromagnetic exchange interaction and two kinds of the 
ferromagnetic interactions; One is induced by the Na- 
gaoka mechanism, which leads to long-range ferromag- 
netic fluctuations for slightly-doped systemS ) 33 i 34 where 
the mechanism is known to work even for finite dop- 
ing levels! 36 ' 55 The other is the three-site ring-exchange 
mechanism, which yields ferromagnetic spin correlations 
for the intermediate filling^ This mechanism works only 



when the product of three hopping integrals along the 
triangles forming the triangular lattice is positive, i.e., 
t\t2 > in our system [see, Fig. [2(a)]. Away from half 
filling, the spin gap has so far been calculated with the 
DMRG methods for some parameters ) 16 ' 19 ' 54 and we now 
study the spin gap in a wider range of the n—t^jtx plane. 
The results are the following: 

(i) For t%/ti = —1, n < 1 [Fig. Efa)], the spin gap is 
considerably enhanced with decreasing n near half fill- 
ing. Since the geometrical spin frustration is reduced by 
doping, the spin-singlet bound state is in a better po- 
sition to be formed. The value of A s increases rapidly 
as 11 decreases, reaches the maximum value A s /ti ~ 0.2 
around n ~ 0.8, and goes down to zero at the FP bound- 
ary n ~ 0.64. The spin gap is always zero when the 
system has two Fermi points, which is in agreement with 
the weak-coupling phase diagram. We should note that 
no singularity in the spin gap is found at half filling. 

(ii) For t 2 /ti — —2, n < 1 [Fig. EJc)] , the spin gap is 
enhanced by doping in the vicinity of half filling as in 
the case of (i). However, unlike in the case (i), the spin 
gap vanishes around n ~ 0.93 even though the system 
has still four Fermi points. This may be related to the 
ferromagnetic spin fluctuations induced by the Nagaoka 
mechanism. Thus, the region for 0.81 < n < 0.93 is 
spin gapless. Further away from half filling, the Nagaoka 
mechanism can no longer work well and the ferromag- 
netic fluctuations weaken. Consequently, the spin gap 
opens again in the filling from n ~ 0.81 to the FP criti- 
cal boundary n ~ 0.45. It is also interesting to note that 
the critical filling n ~ 0.81 agrees with the TLL critical 
boundary; i.e., the spin gap starts to open at the point 
of K p = 0.5. 

(iii) For £2/^1 = 1, n < 1 [Fig- H{b)], the spin gap is 
always finite for small U values [see Ref|16(. For large 
U values, however, the spin gap behaves intricately as a 
function of doping due to the existence of two types of the 
ferromagnetic fluctuations. Near half filling, A s decreases 
with decreasing n and disappears around n ~ 0.88. We 
can see that no spin gap exists in a region 0.57 < n < 
0.88. This is consistent with the fact that the critical 
interaction strength U c at the ferromagnetic transition is 
relatively small for this region^ which would indicate 
the strong ferromagnetic fluctuations there. And then, 
with decreasing n, the spin gap opens again around n ~ 
0.57, where the Nagaoka mechanism no longer works well. 
Like in the case (ii), the point where the spin gap opens 
is on the TLL critical boundary. By further doping, the 
spin gap closes around n ~ 0.4. For n < 0.4, the spin gap 
is zero due to the ferromagnetic fluctuations induced by 
the three-site ring-exchange interaction. Note that there 
is no spin gapless region derived by the three-site ring- 
exchange interaction in the case of (ii) where t\t2 < 0. 

(iv) For t 2 /ti = 2, n < 1 [Fig. H^d)] , the properties are 
qualitatively the same as in the case of (iii). The spin 
gap remains finite only for a tiny region in the vicinity 
of half filling (n > 0.97) and for a small region adjacent 
to the TLL critical boundary (0.74 < n < 0.85). The 
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FIG. 7: Calculated values of the total-spin quantum number 
S as a function of t2/£i for various system sizes L. U/ti = 10 
and n — 1/6 are assumed. 



three-site ring-exchange interaction is more robust than 
in the case of (ih)r^ so that the gapless area becomes 
wider (n < 0.74). 
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E. Fully polarized state 

Of further interest is the presence of the fully polarized 
ferromagnetic state, which occurs when t% is positive in 
the strong-coupling regime. Previously, for U — oo, fer- 
romagnetism has been analytically shown to exist in the 
three limiting cases: n — ► i 2 — > an d n 0i2£. 
Also, for finite U values, it has been shown numerically 
that there is an extensive ferromagnetic phase, where 
the exact diagonalization^ variational^ and DMRG — 
methods have been used. 

Let us then investigate how the ferromagnetic phase 
appears in the phase diagram. We can find it by calcu- 
lating the expectation value of total-spin operator S in 
the ground state, which is defined by 

= £ =S(S + l). (25) 

ij 

For a fully-polarized state, one will obtain S — S'max = 
N/2, i.e., S/ S'max = 1- In Fig. we show the total spin 
S normalized with respect to 5 max as a function of ti jt\ 
at U = 10 and n = 1/6 for various system sizes. We can 
see a transition between paramagnetic and ferromagnetic 
states at (iaAi)c ~ 0.1 and 0.95. The change in S/S max 
at (<2/*i)c becomes sharper with increasing L, suggesting 
the transition to be (almost) discontinuous for L > 72. 
Therefore, the transition may be of the first-order in the 
thermodynamic limit. Thus, the critical transition point 
can be determined in the parameter space, which will be 
given in the next subsection. 



FIG. 8: (a) Ground-state phase diagram of the ID Hubbard 
model with the next-nearest-neighbor hopping, calculated by 
the DMRG method. U/ti = 10 is assumed. The bold (dotted) 
line indicates the FP (TLL) critical boundary and the shad- 
owed area represents a spin-gapped phase. We use the follow- 
ing abbreviations; PM: paramagnetic metal, SG: spin-gapped 
liquid with incommensurate spiral correlations, f-PM: param- 
agnetic metal with strong ferromagnetic fluctuations, FPFM: 
fully polarized ferromagnetic metal, and SS (TS): spin-singlet 
(triplet) superconductivity, (b) Boundary lines between the 
TLL and non-TLL-like regions. 



F. Phase diagram 

Based on the calculated results of the TLL parame- 
ters K p r7 , spin gap A s , and total-spin quantum number 
S, we draw a U = 10 phase diagram of the ID Hub- 
bard model with the next-nearest-neighbor hopping. The 
result is shown in Fig. |5Ja). We find that our system 
Eq. (TT|) exhibits a variety of phases in the parameter 
space of t%/ti and n; it includes a paramagnetic metal 
(PM), a spin-gapped liquid with incommensurate spiral 
correlations (SG), a paramagnetic metal with strong fer- 
romagnetic fluctuations (f-PM), a fully-polarized ferro- 
magnetic metal (FPFM), a spin-singlet superconductiv- 
ity (SS), and a spin-triplet superconductivity (TS). The 
bold and dotted lines in Fig. [E^a) indicate the FP and 
TLL critical boundary, respectively, and the shadowed 
area represents a spin-gapped phase. 

When the system has two Fermi points, the ground- 
state properties are essentially the same as those of the 
standard ID Hubbard model with = 0. The system 
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is thus a paramagnetic metal with 0.5 < K p < 1, where 
the 2/cp-SDW correlation is most dominant. The intro- 
duction of t<z brings a sort of frustration to the 2/cp-SDW 
oscillation, but the oscillation is never broken down as 
long as the system has two Fermi points. 

We now turn to the region with four Fermi points. 
As shown in Sec. IV. D, this region consists of the TLL 
{K p > 1/2) and non-TLL-like (K p < 1/2) phases [see 
Fig.HlJb)]. Looking first at the TLL phase with K p > 1/2, 
the superconducting correlation is most dominant ac- 
cording to the TLL theory. The superconducting phase 
for t2 > is further divided into a couple of phases, de- 
pending on the presence or absence of the spin gap; the 
spin-gapless phase extends over a wide range for large 
^2, which is in contrast to the weak-coupling phase dia- 
gram. This spin-gapless area seems to be expanded by 
increasing the on-site interaction U, as compared to the 
previous DMRG results for U = 2tiM- [I n R e f|Io|, this 
phase is characterized as 2 x C1S1 = C2S2.] The ground 
state is featured as the spin-triplet superconductivity, as 
has been confirmed numerically^ On the other hand, 
the spin-gapped phase is characterized by the spin-singlet 
superconductivity, which is remnant of a wide CI SO re- 
gion in the weak-coupling phase diagram. It is partic- 
ularly worth noting that the TLL phase for t-i < is 
always spin-gapped and the spin-singlet superconduct- 
ing correlation is most dominant. This is consistent with 
the fact that the three-site ring-exchange interaction for 
spin-triplet coupling does not work if t\t2 < 0. 

The other TLL phase belongs to the fully-polarized 
ferromagnetic metal near the FP boundary and at low 
densities. A nearly flat-band system is realized since the 
two band minima arc slightly occupied by electrons at low 
densities (or the band maximum at k = is slightly occu- 
pied by holes near the FP boundary). Consequently, the 
ferromagnetic ground state is stabilized. In the FPFM 
phase, we estimate the TLL parameter as K p ~ 0.5 [see 
Fig. Efb)] , which is the same as that of a spinless fermion 
system. 

We next consider the non-TLL-like regime, which ex- 
tends between the TLL region and the half-filling line. 
The paramagnetic phase with strong ferromagnetic cor- 
relation is located in the vicinity of the TLL regime, 
where the ferromagnetic fluctuations are enhanced due 
to the Nagaoka mechanism; this phase is denoted as f- 
PM in Fig. [7J At present value of U, the total spin of 
the ground state is zero in the entire area of the f-PM 
phase. As U increases, the f-PM phase is enlarged and 
the system would be fully polarized at a critical value of 
U. 

Further approaching the vicinity of n = 1, the spin- 
gapped region appears again. Although most of the spin- 
gapped region is paramagnetic, a narrow spin-singlet su- 
perconducting phase with K p > 0.5 exists along the FP 
boundary line. We can interpret this situation by assum- 
ing the system to be a slightly doped J1-J2 Heisenberg 
model (|24[) : for J2/J1 ^ 0.5, the ground state is of the 
Majumdar-Ghosh state with incommensurate spiral cor- 



relations and the spin-singlet bound state is formed along 
the t2-chains, where the spin-singlet bound state cannot 
easily move. On the other hand, for 0.241 < J2/J1 0.5, 
the ground state is of a dimerized zigzag-bond state 
where the spin-singlet bound state is formed between 
the two ^-chains. At finite doping levels of holes, the 
spin-singlet pairs are mobile in this region, so that in the 
ground state an additional pair of holes is actually con- 
fined to a 'rung' because the gain in kinetic energy due 
to the hole motion is larger than the combined loss in 
the pairing energy and kinetic energy of the spin dimers 
in the Majumdar-Ghosh state. Thus, the narrow spin- 
singlet superconducting state can be regarded as the 
doped zigzag-bond state. 



V. SUMMARY 

We have studied the ID Hubbard model with the 
nearest-neighbor and next-nearest-neighbor hopping in- 
tegrals by using the DMRG method and HF approxima- 
tion. Based on the calculated results of the TLL pa- 
rameters, spin gap, and total-spin quantum number, we 
have determined the ground-state phase diagrams in the 
weak-coupling (U — 2ti) and strong-coupling (U = 10ii) 
regimes. Surprisingly, the strong-coupling phase diagram 
contains a large variety of distinct phases, depending on 
the hopping integrals and band filling. 

We have found for U = 2t\ that the HF results agrees 
well with the DMRG results except for ti as and n as 1 
where the umklapp scattering strength is dominant. The 
phase diagram is qualitatively the same as that in the 
weak-coupling limit. When the system has two Fermi 
points, the 2fcF-SDW correlation is most dominant with 
1/2 < K p < 1. As soon as the Fermi surface splits from 
two into four points, the parameter K p drops (almost) 
discontinuously to 1/2. We then have found K p > 1/2 
in the entire region with four Fermi points and thus the 
superconducting correlation is most dominant. 

We have then found for U — 10ii that the HF re- 
sults no longer agree with the DMRG results because the 
umklapp scattering strength becomes very large. Due 
to the unconventional combination of the charge and 
spin degrees of freedom induced by the next-nearest- 
neighbor hopping, the system can accommodate a va- 
riety of physical phenomena unparalleled in the simple 
ID Hubbard model. The region with two Fermi points 
is characterized by the 2fcp-SDW phase as in the case of 
U = 2t\. However, the region with four Fermi points 
is drastically affected by the Coulomb interaction and 
the breakdown of the TLL state occurs near half filling. 
The strong-coupling phase diagram contains a large num- 
ber of distinct metallic phases; namely, a paramagnetic 
metal, a spin-gapped liquid with incommensurate spiral 
correlations, a paramagnetic metal with strong ferromag- 
netic fluctuations, a fully-polarized ferromagnetic metal, 
a spin-singlet superconductivity, and a spin-triplet super- 
conductivity. 
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In contrast to the weak-coupling phase diagram which 
predicts a spin-gapped liquid (or superconducting) phase 
when the system has four Fermi points, we have found 
that the spin gap vanishes in the substantial region in 
the strong-coupling phase diagram. The absence of the 
spin gap is derived by three types of ferromagnetic mech- 
anisms. The first is the flat-band mechanism around the 
FP boundary and at low densities. A nearly flat band is 
realized since the two band minima are slightly occupied 
by electrons at low densities (or the band maximum at 
k = is slightly occupied by holes near the FP bound- 
ary.) Thus, the ground state is stabilized as the fully 
polarized ferromagnetic state. The second is the three- 
site ring-exchange interaction for i\ at intermedi- 
ate filling, where all the triangles formed by the hopping 
integrals satisfy the ferromagnetic sign rule t\ti > 0. 19 
The ferromagnetic interaction is short ranged^ and is 
not sufficient to make the system ferromagnetic. The 
third is the Nagaoka mechanism near half filling. Since 
the competing antiferromagnetic exchange interaction is 
large, the ground state is not fully polarized. 

Concerning the superconductivity, we have found a 



couple of new phases, which are absent in the weak- 
coupling limit. One is the spin-triplet superconducting 
phase for £2 > ti, where the attractive interaction is 
caused by the gain in kinetic energy due to ring exchange 
of electrons. The other is the spin-singlet superconduct- 
ing phase along the FP boundary near half filling, where 
the superconducting fluctuations are enhanced by large 
difference between two Fermi velocities. This state is also 
regarded as the doped zigzag-bond spin-gapped state. 
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